A superformula for neutrinoless double beta 
decay II: The short range part 



H. Pas ^, M. Hirsch ^, H.V. Klapdor-Kleingrothaus ^, 

S. G. Kovalenko ^ 

^ Max-Planck-Institut fiir Kernphysik, 
P.O.Box 10 39 80, D~69029 Heidelberg, Germany 

^ Inst, de Fisica Corpuscular - C.S.I.C. - Dept. de Fisica Teorica, Univ. de 
Valencia, 46100 Burjassot, Valencia, Spain 

^ Departamento de Fisica, Universidad Tecnica Federico Santa Maria, Casilla 

110-V, Valparaiso, Chile 



Abstract 

A general Lorentz-invariant parameterization for the short-range part of the Qu(3f5 
decay rate is derived. Combined with the long range part already published this 
general parameterization in terms of effective B — L violating couplings allows one 
to extract the Qvl3f5 limits on arbitrary lepton number violating theories. 
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1 Introduction 



The search for neutrinoless double beta decay has been proven to be one 
of the most powerful tools to constrain B — L violating physics beyond the 
standard model [1]. In recent years, besides the most restrictive limit on the 
effective neutrino Majorana mass [1-3], stringent constraints on several the- 
ories beyond the Standard Model such as i?-parity violating [4-6] as well as 
conserving [7] SUSY, leptoquarks [8] and left-right symmetric models [9] have 
been derived (for a review see [1]). While the neutrino mass limit is based on 
the well-known mechanism exchanging a massive Majorana neutrino between 
two standard model V — A vertices, the effective vertices appearing in the new 
contributions involve non-standard currents such as scalar, pseudoscalar and 
tensor currents. Moreover, besides contributions with a light neutrino being 
exchanged between two separated vertices, the so-called long-range part, a.d- 
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ditional contributions from the short range part are expected, e.g. in SUSY 
without R-parity [10,11]. 

Thus we felt motivated to consider the neutrinoless double beta decay rate in 
a general framework, parameterizing the new physics contributions in terms 
of all effective low-energy currents allowed by Lorentz-invariance (see Fig. 
1). Such an ansatz allows one to separate the nuclear physics part of double 
beta decay from the underlying particle physics model, and to derive limits on 
arbitrary lepton number violating theories. The first step of this work, treating 
the long-range part (contributions a-c in Fig. 1), has been pubhshed in [12]. In 
the following the remaining part, treating the short range part (contribution 
d), is presented. 

Although the general decay rate is independent of the underlying nuclear 
physics model, to extract quantitative limits, values of nuclear matrix elements 
are needed, which have been calculated in [11] in the (pn)QRPA aproach. 
Using these, for the first time limits for all potential new physics contributions 
are derived here. 



2 General Pcirametrization 

In the short range part the effective interaction can be considered as point-like, 
thus the decay rate results from first order perturbation theory, 

Roupp^i J d^x{{Z + 2,A),2e-\Cix)\iZ,A),-). (1) 

The most general Lorentz invariant Lagrangian has the following form 

+eeJ^rj^. + e^jr^j^, + e^J^^r'^j'^}, (2) 

with the proton mass rUp and the hadronic currents of defined chirality J = 
u{l ± 75)d, J^' = uYil ± lb)d, J^" = m^f^^^YK^ ± 75)c^ and the leptonic 
curents j = e(l ± 75)6'^, = e7''(l ± 75)6'^ , jf^" = e^[^>^, YK'^ ± Ts)^*^- In 
some of the cases the decay rate for the effective coupling may depend also 
on the chirality of the currents involved. In these cases we define = e^^^, 
where xyz = L/ R, L/ R, L/ R defines the chirality of the hadronic and leptonic 
currents in the order of appearance in eq. (2). In the cases where it is not 
necessary to distinguish the different chiralities we suppress this additional 
index. 

In renormalizable theories no fundamental tensors exist. Thus the tensor cur- 
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rents have to result either from Ficrz rearrangements or by integrating out 
heavy particles, when deriving the effective Lagrangian from the fundamental 
theory and decomposing the expressions obtained in terms of the Lorentz in- 
variant bilinears used above, e.g. from ■U7^7,^( 1 + 75)0? = g^yJ — iJ^u- Moreover, 
the contribution of the leptonic tensor current vanishes in s-wave approxi- 
mation. Since the tensor currents in the rearrangements and decompositions 
mentioned above are always created together with dominant contril:)utions, the 
corresponding terms proportional to eg, €7, eg can be neglected. The remaining 
terms are the ones proportional to ei, £2, €3, €4, £5. 

Applying the standard nuclear theory methods based on the non-relativistic 
impulse approximation we derive the general Oi//3/3-decay half-life formula in 
s-wave approximation 

K^'r - (3) 

Gil Ee^Mr + G2\j2^^M^\' + GsRe[if2^Mi)ijZeMin 

i=l i=4 i=l i=4 

Here the phase space factors are 

Gi — Goi, G2 — - — ^-^Go9, G3 = ^ ^ ^ G06, (4) 

with Gok given in [13] and shown in Table 1. 
The nuclear matrix elements in eq. (3) are 

Ml = -aiMF,N, M2 = -a2MGT,N, 

Ms = —{Mgt,n T asMF,N}^ (5) 

mpme 

Mi = ±aiMGT,N, M5 = T(^5Mf,N- 

In the last three cases contractions of hadronic currents with different chirali- 
ties lead to different results. The negative sign in M3 corresponds to e^^^ and 
^RRz ( Jy-|-^Jy-|-^), the positive sign to 63-^^ and ef^^ {JvtA'Jv±a)- The sign of 
is positive for the combinations £4^^, ef^^, ef^^, e^^^ {JvtaJtl/trJv-a 
and Jv±AJTL/TRjv+A)- For the combinations e^^^, ef^^, ef^^, e^^^ {Jv^aJtl/trJv+a 
and Jv±AJTL/TRjv-A) it is negative. The sign of the matrix element M5 is 
negative for the left-handed leptonic current (ef^^) and positive for the right- 
handed one {e^^^). 

The standard nuclear matrix elements M f,n und M.gt,n in eq. (5) are defined 
as 



Mgt,n = (F| E^Vi^V, • a, (^] F^ixAm, 



(6) 
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Mp,N = {F\^rM^ ( ^ ) FivMI/). (7) 

Here Vij and Ro are the distance of the nucleons involved and the nuclear 
radius respectively. Numerical values of these matrix elements for '''^Ge [11] 
in the Quasi Particle Random Phase Approximation(pn-QRPA) are given in 
Table 1. Uncertainties of nuclear matrix element calculations are discussed in 
Ref. [11,14]. 

The prefactors ctj in eq. (5) are defined as follows, 



QA ' mprUe ^ QA ' mpme 



f9v^^ m\ gyF^ 



as = I , 0:4 = Q;5 = „ 

gA mpuie Qa mpme 

The finite nucleon size is taken into account in a common way [15] by intro- 
ducing the nucleon form factors in a dipole form 

9v,A Fs Tf [ ' 

with niA = 0.85 GeV, gy = 1-0, (7a = 1-26. The other form factor normaliza- 
tions have been calculated in Ref. [16] within the MIT bag model, 

fP = 0.48, t/^' = 1.38. (10) 

The momentum dependence has been absorbed into the integral over the mo- 
mentum q transferred between two nucleons, 

F^ixA) = 4nmln,J ^^^y ^^2^\ ^y ^'""'' = ^(3 + 3xA + x^)e--,(ll) 
where xa = mATij- 



3 General 0i//9/9-decay constraints 



Let us extract limits on the Icpton-number violating parameters in eq. (3) 
from the experimental lower bound on the 0z//9/9-decay half-life. The conser- 
vative half-life limit of the Heidelberg-Moscow experiment [2] 

T°f ^ > 1.8 • 1025years(90%aL.) (12) 
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provides at present the most stringent bound. Applying this hmit to eq. (3) 
leads to a complex 5-dimensional exclusion plot for all five parameters e^. To 
make the Oi//3/3-decay experimental bound more transparent it is a common 
practice to assume the dominance of one coupling. This means only one non- 
zero Ei is considered at a time (evaluation "on axis"), while the interference 
between different contributions is neglected. In this way we obtain the upper 
bounds on listed in Table 2. 

In the general case, when more than one Ea is involved simultaneously, we deal 
with the multidimensional exclusion plot. For example, this situation occurs in 
the left-right symmetric models, see [9,13]. It also appears natural, if effective 
couplings are generated in Fierz rearrangements or decompositions of given 
expressions. In the following, as an example, we dicuss the case of interfer- 
ence of e^^^ and e^^^ contributions (the cut in the ef^^-Cg^'^-plane), which 
can result from the decomposition '^7^7,^(1 -|- ^b)dJ'^j^ — {g^uJ — iJ^iv)J'^3^ 
mentioned above. In this case the expression for the inverse half life becomes 

[Tfi^^V = \e'i''''M, + e',^^M,\'G,, (13) 

leading to the exclusion plot shown in Fig. 2. The region outside the ellipse 
is excluded. It is obvious that the limits on e^^^ and 65^"^ in this case are 
somewhat less stringent (by about 15-20%) than the on-axis limits listed in 
Table 2. The analogous treatment of other combinations, using the matrix 
elements hsted in Table 1, is straightforward. 



4 Conclusion 

We have presented a general parameterization for the short range part of 
the neutrinoless double beta decay rate in terms of effective couplings. The 
resulting bounds are summarized in Table 2. Combined with the long range 
part already published [12], this parameterization provides the double beta 
decay constraints for arbitrary lepton number violating theories beyond the 
SM. 
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Fig. 1. Feynman graphs of the general double beta rate: The contributions a) - c) 
correspond to the long range part discussed in [12], the contribution d) is the short 
range part. 
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Fig. 2. Simultaneous treatment of and . The Oi^/3/3 half life limit allows 

for a region (inside the ellipse) in a two dimensional parameter space. The limits 
on e^^^ and e^^^ in this case are somewhat (by about 15-20%) less stringent than 
the ones evaluated "on axis". 
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Mgt,n 


Mf,n 


Goi(yr-)"' 


Go6(yr)-i 




1.13 • 10-1 


4.07 • 10-2 


6.4 • 10-15 


1.43 • 10-12 


3.3 • 10-1° 



Table 1 

Nuclear matrix elements for ^^Ge 0u(5(5-decay calculated in the pn-QRPA approach 
(from [11]) and phase space factors (from [13]). 





k2| 




\^LRz\ \^RLz\ 
l'F3 1 


k4| 


k5 


3 • 10-^ 


2 • 10-^ 


4 • 10-s 


1 • 10-8 


2 • 10-8 


2 • lO-'^ 



Table 2 

Experimental upper bounds on the absolute values of the effective B — L violating 
couplings evaluated "on axis". For 63 the contractions of hadronic currents with 
different chiralities lead to different results. 
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